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Abstract. It is shown that the layered antiferromagnetic order in stoechiometric LaMnOs cannot be un- 
derstood purely from electronic interactions. On the contrary, it mainly results from strong cooperative 
Jahn-Teller deformation. Those involve a compression of the Mn — O octahedron along the c— axis (mode 
Qz < 0), while alternate Jahn-Teller deformations occur in the ab— plane (mode Qz)- These deformations 
stabilize a certain type of orbital ordering. The resulting superexchange couplings are calculated by exact 
diagonalization, taking into account both e g and t2 g orbitals. The main result is that antiferromagnetic 
(ferromagnetic) coupling along the c— direction (ab— planes) can be understood only if the Jahn-Teller 
energy is much larger than the superexchange couplings, which is consistent with experiments. This mech- 
anism contrasts with that based on weak Jahn-Teller coupling which instead predicts elongation along the 
c— axis (Q3 > 0). The crucial role of the deformation anisotropy ^ is also emphasized. 

PACS. 71-70. E Jahn-Teller effect - 75-30. ET Exchange and superexchange interactions 



1 Introduction 

Perovskite oxides containing Mn ions have been the ob- 
ject of intense interest in the recent years. In spite of being 
known for a very long time, these compounds have been 
reconsidered in great detail owing to their colossal magne- 
toresistive properties. Starting from the "parent" phases 
LaMn0 3 (trivalent Mn) and CaMn0 3 (tetravalent Mn), 
substitutional doping has revealed an extremely rich phase 
diagram. Understanding this diagram requires at least the 
following ingredients : i) strong on-site Coulomb interac- 
tions; ii) the "double exchange" mechanism due to the in- 
terplay of e g electron itineracy and Hund's exchange with 
the more localized t 2g electron spins, which favours fer- 
romagnetism [fil0||; iii) superexchange between t 2g elec- 
trons as well as between e g electrons on neighbouring sites; 
iv) large electron-lattice interactions, inparticular due to 
Jahn-Teller (JT) effect on Mn 3+ ions j§|]. All these ele- 
ments are necessary to understand the interplay between 
spin, charge and orbital ordering. The latter lifts the de- 
generacy of the e g orbitals by a cooperative Jahn-Teller 
lattice deformation and leads to tetragonal or orthorhom- 
bic deformations of the cubic structure. 

Although Goodenough || provided long time ago a 
qualitative understanding of the phase diagram of the 
(La,Ca)Mn0 3 family, a full microscopic description is 



still lacking. Especially the dramatic dependence of all 
physical properties with very fine tuning of the chemi- 
cal composition requires a precise estimate of the various 
parameters, and clear identification of the dominant mech- 
anism for every doping. Surprisingly enough, such an un- 
derstanding is not yet reached in the insulating antiferro- 
magnet LaMnO^, although it seems essential before quan- 
titatively studying the doped phases. This phase, when 
fully stoechiometric, presents a layered antiferromagnetic 
order, with ferromagnetic couplings (F) in two directions 
and antiferromagnetic ( AF) coupling in the other . The 
AF directions are associated to a shortening of the Mn—O 
bonds, leading to tetragonal distortion, while in the F di- 
rections long and short bonds alternate, yielding the over- 
all orthorhombic structure. In what follows, we shall ne- 
glect the tilting of the Mn — O octahedra and concentrate 
only on the Mn — O bond length deformations. These can 
be understood in terms of cooperative JT effect. The cor- 
responding lifting of e g degeneracy can be viewed as an 
orbital ordering, with occupied d orbitals pointing prefer- 
entially in the directions of long Mn — O bonds. 

Several proposals have been made to explain layered 
antiferromagnetism in LaMnOs . Goodenough || used the 
picture of " semi-covalence" where oxygen orbitals play an 
essential role in overlapping empty d orbitals of Mn ions. 
This picture, although useful for qualitative purposes, has 
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not received confirmation by microscopic calculations and 
does not allow to write simple enough models, for instance 
based on a Hamiltonian involving only metal orbital elec- 
trons and their basic interactions. A microscopic descrip- 
tion requires to identify clearly the dominant interactions 
in the problem. In pioneering works, Kugel and Khom- 
skii B, and Lacroix || (see also earlier work by Roth 
ifTof), proposed that superexchange in the presence of e g 
orbital degeneracy results in ferromagnetism and orbital 
ordering : Hund's rule favours in this case different or- 
bitals on neighbouring sites and ferromagnetic coupling. 
Using a simplified model with equal hopping integrals be- 
tween e g orbitals leads to the same ordering along the 
three cubic lattice directions : the resulting structure is 
an insulating ferromagnet, with " antiferroorbital" order- 
ing. However, taking properly into account the hopping 
integrals between d x 2_ y 2 (denoted x) and d z 2 (denoted z) 
orbitals, Kugel and Khomskii j|] found the correct mag- 
netic structure. Starting with degenerate e g orbitals, they 
performed a perturbative calculation in jj and ^ff- where t, 
Jh and U are the typical hopping integral, the Hund cou- 
pling and the on-site repulsion in the order. Based on the 
weak electron- lattice coupling in the compound KCuF^, 
they considered the JT couplings as a perturbation. As 
a result, orbital and magnetic ordering result from su- 
perexchange (SE) only: Intraorbital SE dominates in the 
c-direction (defined as the z-axis), leading to AF coupling, 
while interorbital SE dominates in the ab-directions, yield- 
ing F coupling. Occupied orbitals are dominantly d z 2_ x 2 
and d z 2_ y 2 1 therefore, as Kugel and Khomskii remark, 
for Cu 2+ in KC11F3 (hole orbital), JT coupling implies 
a shortening as experimentally observed (c/a < 1). How- 
ever, for Mn 3+ ions with one electron in the e g levels, they 
correctly point out that repulsion between metal and an- 
ion orbitals, together with JT coupling, would trigger a 
lengthening of the c-axis (c/a > 1), in contradiction with 
the actual structure. In a recent work, Feiner and Oles 
[ p"l[ reconsidered Kugel and Khomskii's model, including 
both Hund's coupling betwen e g and t 2g orbitals and the 
antiferromagnetic superexchange interaction between t 2g 
spins (equal to § in the ground state). Their results con- 
firms those of Ref . || : They find the correct layered struc- 
ture (which they call MOFFA) , but only if the d z 2 orbital 
has lower energy than the d x 2_ y 2 one, contrarily to what 
happens for electron-like orbitals (case of LaMnO^). This 
contradiction sets the limits of the Kugel-Khomskii model 
for LaMnOz. We believe that the JT effect, on the con- 
trary, has to be considered from the very beginning in the 
model. 

Essentially, the assumption that the e g degeneracy is 
lifted principally by superexchange may be justified in 
KC11F3, but is definitely not correct in LaMnO^. In fact, 

this could hold only if the typical JT splitting e was much 

.2 

smaller than the superexchange splitting, of order jj . The 
latter (related to the magnetic transition temperatures) 
being of the order of a few meV , the former is much larger. 
Although there is no precise evaluation of this quantity, 
this is supported by experiment: On the one hand, the de- 
formations of Mn-0 bonds is extremely large, more than 



ten per cent, indicating that e > ksT. On the other hand, 
neutron scattering measurements show that the local dis- 
torsions per sist above the orthorhombic-cubic transition 
at 750 A" |12f . This temperature only marks the disappear- 
ance of cooperative JT ordering, while distorted MuOq 
octahedra still exist at higher temperatures. Photoemis- 
sion [T^ | measurements indicate that JT splittings are as 
large as a few tenth of eV, comparable to the electronic 
hopping integrals between neighbouring sites. And opti- 
cal conductivity analysis [jl4| also shows evidence of large 
splittings. 

In these conditions, the degenerate perturbation calcu- 
lation of Ref. |^| does not hold anymore. In a previous work 
we have reconsidered the problem wihin perturbation the- 
ory pB| , making the opposite assumption, i.e. e >> |j: 
This means that, given the crystal deformations, due to 
strong cooperative JT effect, the e g orbitals split so as to 
give a certain type of orbital ordering. The orbitals sta- 
bilized at each sites are different from the one predicted 
by pure superexchange. We have found that, depending 
on the values of the two JT modes Q 2 and Q3, different 
magnetic ordering could be stabilized, among which the 
layered "FFA". This ordering is always stabilized if the 
Q3 mode is positive, e.g. for dilatation in the c-direction. 
But in the real case Q3 < 0, FFA order is realized only 
if the in-plane alternate Q 2 mode is sufficiently large and 
overcomes the contrary effect of Q3. Looking at structural 
numbers, one checks that this is actually the case. Nev- 
ertheless the system is close to the point where the FFA 
order becomes unstable towards FFF. This results in the 
F exchange (along the ab-plane) being larger than the 
AF exchange (along the c-axis). This feature has been ob- 
tained from inelastic neutron scattering , and it cannot 
be explained by the Kugel-Khomskii model, which obtains 
on the contrary that the F superexchange is of order ^j- 
times the AF one, thus much smaller. 

The interplay between lattice distortions and magnetism 
has also been investigated from ab-initio calculations of 
the electronic structures |l7| , p^|Jig[| . All conclude with a 
prominent role of those distortions to stabilize the actual 
magnetic order. In particular, Solovyev et al. have 
found that the c-axis exchange is antiferromagnetic only 
if the JT distortion is sufficiently large. For LaMnOz with 
its very large distortion they obtain the layered antiferro- 
magnetic structure, but it is close to the border between 
FFA and FFF phases. Very recent Monte Carlo calcu- 
lations have also demonstrated the relevance of the JT 
interaction in stabilizing the FFA magnetic order [^o| . 

In the present work, we reconsider the problem, be- 
yond any perturbation theory, by exact diagonalizations 
on pairs of Mn 3+ sites. The two e g orbitals are considered 
together with the quantum |-spins due to the electrons in 
the t 2g levels. Our conclusions confirm the essential role 
of JT deformations, especially the Q 2 mode, to stabilize 
the layered AF order. They also demonstrate that it is 
essential to include Hund's coupling with t 2g orbitals, and 
that the role of the intrinsic t 2g AF exchange is to slightly 
stabilize the FFA order with respect to the FFF one. 
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2 The Model 

^From the discussion of the preceding section it is clear 
that the basic physical ingredients required for a satisfac- 
tory description of the Manganites should involve both 
Coulomb and lattice (namely JT) interactions. Accord- 
ingly we consider the following model 

H = H t + H H + Huu, + Hj + H JT (1) 

with 

iaaa' er 

Hh = —Jh c icea s aa tCi a<T f 

a' ' ^oloct' 

Huu> = U J2 ( c L T c «t) ( c Li c «*l) 

iot 

+ U ' ( C La c ^a) (cLv'CiaV') 

ia^a' aa' 

Hj = J t ■ Sj 

(«•> 

H.JT = 9^ ( C La r(3 ' ] aa'Cia'aQzi + c\ aa T (2) aa : C la > <T Q 2 ij 

i 

The first term represents the kinetic energy with the 
electrons in the Manganese 3d x 2_ y 2 (a = x) or 3d3 2 2_ r 2 
(a = z) orbitals hopping from site i to the nearest neigh- 
bor (nn) site i + a in the a lattice direction. Here s is 
the vector of Pauli matrices for spins and r the vector of 
Pauli matrices for orbital pseudospins in the x, z basis. 
Specifically, for a standard choice of the phases for the or- 
bital wavefunctions, the hopping between the x and the z 
orbitals are given by 

^x'x 3i, t z '^ — t; 
t* z = -V3t; t y xz = V3t 
t zz = -U t% x =tl z = 0. (2) 

Together with the Hund coupling given by Hjj the ki- 
netic energy gives rise to the usual "double-exchange" itin- 
erancy of the e g electrons. The (strong) on-site Coulomb 
interactions, are represented by the intraorbital repulsion 
U and by the interorbital U' = U — 2Jh term. 

For simplicity here and in the following we will not 
distinguish between the Hund exchange energy between 
electrons in the e g and ti g orbitals. The antiferromagnetic 
superexchange coupling between neighboring ti g spins is 
considered with Hj, while the JT interaction between the 
e g electrons and the (cooperative) lattice deformation is 
given by the last term Hjt- The Jahn- Teller modes are 
defined in terms of the short (s), medium (m) and long (I) 
Mn-0 bonds by Q 2 = y/2(l - s) and Q 3 = v / 2/3(2m - 



I — s), the m bonds lying in the z direction and the s,l 
ones in the x, y planes. 

Since in the present work we will not attempt to per- 
form any energy minimization by including the elastic in- 
teractions due to the lattice, we here disregard these en- 
ergy terms by treating the JT deformations Q = (Q2, Q3) 
as external fields imposed by a lattice ordering involving a 
much higher energy scale than the magnetic ones. There- 
fore in the following the various magnetic couplings will 
be determined in terms of assigned lattice deformations. 
This viewpoint, which already guided us in the perturba- 
tive analysis of the stability of FFA antifcrromagnetism in 
the undoped LMO [jl5| is definitely justified by the exper- 
imental observation that the JT energy splitting is much 
larger than all magnetic couplings. 

We exactly diagonalize the Hamiltonian in Eq. (|l|) for 
a system of two sites with open boundary conditions. The 
two sites are located cither on the same xy plane or on 
adjacent planes and the suitable hopping matrix elements 
between the various orbitals have been considered accord- 
ing to expressions (||). 

The JT energy splitting e = g\jQ\ + Q\ and the defor- 
mation anisotropy ratio r = Q2/Q3 are given external pa- 
rameters and are fixed for any diagonalization procedure. 
Once the ground state is found, the effective exchange 
coupling between the total spins on the two sites can be 
determined. Specifically, since the Hamiltonian conserves 
the total spin of the two-site cluster, we determine the 
ground states with total spin St = 4, Ms T — 4 and 
St — 3, Ms T = 3. Then the magnetic coupling is given 
by the energy difference E(S T = 4, M St = 4) - E(S T = 
3, Ms T — 3) = 2 J. Once the magnetic couplings (and par- 
ticularly their sign) along the various lattice directions are 
found, the resulting magnetic phase is also determined. 



3 Results 

In order to gain insight from the physical processes un- 
derlying the intersite magnetic couplings, we first carry 
out a comparison between the results of the perturbative 
analysis of the superexchange interactions (see Ref. p5|) 
and the exact numerical calculations. The perturbative 
analysis not only was performed assuming very large local 
Coulomb interactions (U, U' and Jh much larger than t), 
but the additional assumption was made that the JT en- 
ergy splitting e greatly exceeds the typical superexchange 
energy scale of order t 2 /U. In this way the ground state 
can safely be assumed to be formed by just one singly 
occupied e g level. Accordingly the exact numerical cal- 
culations to be compared with the analytic results have 
been performed for Q3 < and t — 0.2eV, U — 8eV, 
J H = 1.2eV, e = OAeV, and J t = 0. Fig. 1 reports the 
superexchange interactions both in the planar and inter- 
planar directions obtained with both the perturbative and 
the exact-diagonalization analysis. As it is apparent, the 
perturbative J xy and J z display the same qualitative be- 
havior as in the exact calculation. This confirms that, at 
least in the e t 2 /U limit, a substantial part of the mag- 
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Fig. 1. Magnetic couplings for Q 3 < 0, t = 0.2eV, U = 8eV, 
Jh = 1.2eV, e — OAeV, and Jt = as function of the defor- 
mation anisotropy ratio from exact diagonalization and per- 
turbation theory. Dashed line: J xy from exact diagonalization; 
Dot-dashed line: J z from exact diagonalization; Solid line: J xy 
from perturbation theory; Dotted line: J z from perturbation 
theory. Negative (resp. positive) values indicate ferro (resp. 
antiferro) magnetic interactions. 



netic effective interactions is generated by the superex- 
change processes due to the hopping of electrons lying 
in the lower e g level on the same or on different nearest 
neighbor e g orbitals. On the other hand, the quantitative 
comparison indicates that the range of stability for the 
FFA phase (i.e. J xy < and J z > 0) is modified. In fact 
a positive J z together with a negative J xy are obtained in 
the exact calculation on a somewhat larger range of lat- 
tice deformation anisotropics (Q2/IQ3I > 2.5). In order to 
establish a tighter connection between the experimentally 
determined J's and the observed deformations, and to in- 
vestigate the role of the various interactions in the model, 
a more systematic analysis is required. Assuming the JT 
interaction to be relevant in stabilizing the FFA phase, 
we investigate the behavior of the exchange constants J xy 
(denoted "intraplane") and J z (denoted "interplane") in 
terms of e and the deformation ratio r. 

Figs. 2 and 3 report J xy and J z as functions of \r\ for 
the Q3 < case (the one relevant for LMO) at a large 
(e ~ 2t) and at a small (e < O.lt) value of the JT split- 
ting respectively. Different values of the Hund coupling 
Jh are considered. One can first observe from Fig. 2 that 
in the large-e case the increase of the Hund coupling shifts 
downwards both the intraplane and the interplane mag- 
netic couplings. This outcome can be rationalized in terms 
of perturbatively generated superexchange processes pro- 
viding AF effective couplings of the form 
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Fig. 2. Magnetic couplings J xy (solid line) and J z (dotted line) 
vs. the deformation ratio r = Q2/IQ3I for negative values of 
Q 3 , for t = O.UeV, U = GeV, e = 0.3eV, J t = 2.1meV and 
Jh = 1.2eV (Jh = 0.9eV) in the upper (lower) panel. 
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Fig. 3. Magnetic couplings (see Fig. 2) vs. the deformation 
ratio r = Q2/IQ3] for negative values of Q3, for t — O.UeV, 
U = 6cV, e = O.OleV, Jt = 2.lmeV and Jh = 1.2eV (Jh = 
0.9eV) in the upper (lower) panel. 



competing with the generated F interaction 



C 



xy ' z U + e-(5/2)J H ' 



(4) 



U+{3/2)J H U + > 



(3) 



The numerical coefficients A, B, and C stem from the dif- 
ferent hopping matrix elements between the different or- 
bitals in the different directions pa. Specifically, while 
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the ^4's are related to the hopping processes between two 
nearest-neighbour lower- lying e g orbitals, the B and C co- 
efficients are due to hoppings between one low-lying and 
one higher JT-split orbitals (this is why the corresponding 
denominators involve e) . The A, B and C coefficients are 
independent from the Coulomb interactions, which only 
determine the energies of the virtual intermediate states in 
the superexchange processes. The above schematic expres- 
sions clearly show that, when Jh is increased for a fixed 
e ^> t, the F spin configuration becomes more favourable, 
since the F coupling become stronger, while the AF inter- 
action weakens. We remark that purely electronic models 
such as in Refs. ma ke use °f degenerate perturba- 
tion theory. Then the orbital splitting is of order of the 
exchange couplings J and therefore those models become 
invalid if e > J, which is the case in LaMnOj,. More se- 
riously, the orbital order resulting from purely electronic 
interactions is at odds with that obtained from the actual 
Jahn- Teller distortions, showing that those distortions do 
not result from an orbital ordering of electronic origin, but 
are on the contrary the mere source of orbital ordering. 

Another quite generic effect, which can be interpreted 
in terms of perturbatively generated superexchange pro- 
cesses is the tendency of J z to acquire a F (or at least a less 
AF) character at low values of Q2/IQ3I (this can also be 
accompanied by an upturn of J z for |r| tending to zero). 
This occurs because for Q3 < 0, the lowest e g level pro- 
gressively loses its 3d 3z 2_ r 2 component: By schematically 
writing the lower and the upper e g states as \a) oc |cc)+?7|z) 
and 1 6) cx —ri\x) + \z) respectively, 77 vanishes with \r\ — > 0. 
Now, the superexchange along z is driven by the interplane 
hopping, which is only allowed between 3d 3z 2_ r 2 orbitals. 
Furthermore one can see |H| that the ferromagnetic su- 
perexchange arises from \a) — ► \b) hoppings, which are 
of order ?y, while the antiferromagnetic coupling is mostly 
generated from intraorbital \a) — > \a) hopping (the A term 
in Eq. (||). Since this latter is of order rj 2 , it is quite natu- 
ral that in the low-|r| region, as i] decreases, the superex- 
change along z is ferromagnetic and vanishes with r\. This 
ferromagnetic tendency is, however, contrasted (and ac- 
tually overcome in Figs. 2 and 3) by the independent AF 
superexchange Jt between the t2 g spins, which becomes 
relatively more important. Of course, when \r\ increases, 
the rf 1 terms in the hopping become relevant, the intraor- 
bital I a) — > I a) hopping starts to dominate and J z even- 
tually becomes (more) positive (i.e. AF). 

As far as the superexchange along the planes is con- 
cerned, at small \r\ this is instead dominated by the large 
hopping between 3d x 2_ y 2 orbitals, which favor the \a) — > 
I a) hopping and, consequently produces an AF magnetic 
coupling. On the contrary, for large |r|, orbital ordering 
implies that the main superexchange contribution comes 
from hopping between different orbitals, thus favouring 
ferromagnetism 

All the above arguments are obviously only valid as 
long as the conditions for the perturbation theory nearly 
hold. On the other hand, the simple perturbative approach 
between non-degenerate states breaks down when e w 
t 2 /U as in Fig. 3 and the interpretation of the results 
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Fig. 4. Magnetic couplings vs. the JT splitting energy e for 
t = O.UeV, U = 6eV, Jt = 2.1meV for Q2/Q3 = -3 (upper 
panel) and Q2/Q3 = 3 (lower panel). Solid line: J xy for Jh = 
1.2eV; Dotted line: J z for Jh = 1.2eV; Dashed line: J xy for 
Jh = 0.9eV; Dot-dashed line: J xy for Jh = 0.9eV. 



is not so transparent. However, the effect of Jh favoring 
ferromagnetism is still present. 

An important difference between the results in Figs. 2 
and 3 is that the FFA phase is generically obtained in a 
broad range of parameters when t>f. In particular, for 
rather realistic values of Jh ~ 5t ~ leV the deformation 
ratios required to generate negative (i.e. F) couplings in 
the xy planes and positive ones in the z direction are quite 
reasonable |r| ~ 2 — 3. The same does not hold in the case 
of small JT splitting, where J xy and J z have the same sign 
(FFF). Therefore a first result is that a sizable e is needed 
in order to obtain both the FFA phase and reasonable 
lattice distortion ratios Qi/Qz- 

This result is also confirmed by the calculation of J xy 
and J z as a function of e, at a fixed value of the defor- 
mation anisotropy ratio r. Figs. 4 and 5 report the values 
of J xy and J z for r = 3 and r — — 3 respectively. While 
the positive r case is generic for perovskite materials with 
the lattice elongated in the z direction (c/a > 1), the lat- 
ter choice is more pertinent to the case of the undoped 
LMO, where c/a < 1. As already discussed by Kugel and 
Khomskii || for a different model and as confirmed by the 
perturbative analysis of Ref. jl5| , the JT deformation and 
the superexchange interactions cooperate when Q3 > 
like in KCuF 3 so that it is not surprising that for all val- 
ues of Jh the FFA is realized over a much broader range 
of e. On the other hand, for Q3 < 0, Fig. 4 shows that 
the conditions for a FFA phase, J xy < and J z > 0, 
are only realized for a smaller range of e values. In par- 
ticular a sizeable minimum value of e is required to have 
an AF coupling along z, while exceedingly large values of 
e (of order Jh) produce an AF coupling also along the 
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planes. Both the minimum and the maximum values of e 
for obtaining the FFA phase increase upon increasing Jg ■ 
However, the maximum value of e increases more rapidly 
and the overall effect is that, increasing Jh, the available 
range in e to obtain an FFA phase is enlarged. Again the 
behavior displayed in the exact calculations reported in 
Fig. 4 can easily be interpreted in terms of the pertur- 
bative superexchange processes schematically represented 
in Eqs. (||) and (^). First of all these expressions at once 
account for the increasing behavior of the couplings upon 
increasing e: While only the interorbital part of J af (the 
contribution proportional to B) decreases upon increasing 
e, the whole ferromagnetic part in Eq. (^) is suppressed 
when e grows, so that the total coupling, although ferro- 
magnetic at small JT energy splitting, eventually vanishes 
and becomes positive. 

Moreover it turns out that, for \r\ > 2 — 3 the hop- 
pings generate smaller A, B, C coefficients in the z direc- 
tion. This accounts for the more rapid rise of J z when e 
is increased. Finally, along the same line of the discus- 
sion of Fig. 2, one can easily observe that an increasing 
Jh strcnghtens the ferromagnetic component and weakens 
the antiferromagnetic one, thus rationalizing the generic 
tendency of all curves to be shifted downwards when Jh 
grows. 

Besides the above specific findings, the occurrence of 
the various magnetic phases can be cast in a phase dia- 
gram at zero temperature illustrating the stability region 
of these phases in terms of the JT energy splitting and 
the deformation ratio. In the light of their richer com- 
plexity and of the present interest in the Manganites, we 
here consider in greater detail the case of Q3 < of rel- 
evance for the undoped LMO, while the Q3 > case is 
only described in the inset of Fig. 5. Figs. 5 and 6 report 
the phase diagram for two different values of the Hund 
coupling. Both phase diagrams display the same qualita- 
tive features. In particular, at moderate and large values 
of e a Neel AAA phase is found for weak planar distor- 
tions (small r). As seen in the discussion of Fig. 2, in the 
very-small-r region, J xy is naturally positive, while the 
superexchange between e g levels along z, although ferro- 
magnetic, is small so that the direct superexchange be- 
tween tig spins may easily dominate and gives rise to the 
AAA phase (see Figs. 2 and 3). As it can be also be seen 
from Fig. 1, it can be checked that the AAA phase is 
replaced by the so-called C-like antiferromagnetic AAF 
phase in the Jt = case. At small-to-intermediate values 
of e, a progressive increase of |r| drives the system towards 
the phase AAF. In this phase J xy keeps its AF character, 
while the negative superexchange between e g levels along 
z is small, but no longer is overcome by Jt- At larger val- 
ues of e the AAF phase is not present, but the intimate 
nature of the AAA phase changes upon increasing In 
particular while at low |r| the AF along z is determined 
by Jt, at larger |r|, the superexchange between e g levels 
along z is itself AF and therefore the t^g superexchange 
contributes, but it is not strictly necessary to the AF cou- 
pling along z. On the other hand, a further increase of |r| 
promotes a F coupling along the planes and leads to the 
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Fig. 5. Zero-temperature phase diagram for t = 0.14eV, U = 
6eV, Jt = 2.1meV and Jh = 1.2eV '. 
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Fig. 6. Zero-temperature phase diagram for t = 0.14eV, U = 
6eV, J t = 2.1meV and J H = l.OeV. 



A-type antiferromagnetism FFA experimentally observed 
in undoped LMO. 

At small values of the JT splitting, the phase diagram 
is prominently occupied by a FFF phase. In this latter re- 
gard, from the comparison of Figs. 5 and 6, the important 
observation can be done that the FFF phase at low and 
moderate e's is greatly stabilized by the increase of the 
Hund coupling Jh, as previously expected. 

Within the present exact numerical treatment of the 
model in Eq. (|j) it is also possible to attempt at "precise" 
estimates of J xy and of J z . As an example, we report here 
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a realistic sets of parameters (among many others) pro- 
viding the values J xy — —0.83 meV and J z — 0.58 meV 
experimentally observed with inelastic neutron scattering 
pq| . Assuming Qg /IQ3 1 = 3.2, a value largely confirmed 
by many groups |f2], we take t = 0.124eV, U = 5.81eV, 
J H = 1.2eV, J t = 2.1 meV, and e = 0.325eV. The quite 
reasonable values of the model parameters needed to re- 
produce the measured magnetic couplings is an indirect 
test of the validity of the considered model. We empha- 
size that the "anomalous" trend \J X y\ > \Jz\ is correctly 
reproduced, and that our fit is relatively flexible concern- 
ing parameters U, Jh or i, provided e is large enough. 

4 Conclusions 

In this paper we presented the results of calculations based 
on the exact diagonalization of a model aiming to describe 
the stoechiometric LaMnO,3. The model includes strong 
local Coulomb interactions as well as a JT coupling be- 
tween the electrons and the Q2 and Q3 lattice deforma- 
tions. 

Despite the smallness of our cluster, we believe that 
our determination of the magnetic couplings not only is 
qualitatively, but also quantitatively significant. This is so 
because, in the presently considered undoped LMO, the 
coherent charge mobility is negligible due to the large on- 
site Coulomb repulsions and to the substantial JT defor- 
mations. As a consequence the magnetic interactions do 
not arise, e.g., from Fermi surface instabilities or other col- 
lective effects, but are rather determined by short-distance 
(incoherent) processes. 

One first relevant result is that, when the Mn — O 
octahedron is compressed along z, a FFA phase is only 
obtained for a sizable (staggered) Q2 deformation of the 
planar unit cell. This finding agrees with the ab initio 
calculations of Ref. (lq] . 

Our analysis also points out the relevant role played 
by the Hund coupling, which gcncrically emphasizes the 
ferromagnetic component of the superexchange processes. 
Quite relevant turns out to be also the Hund coupling 
between the e g electrons and the ti g spins. In this latter 
regard, we explicitely checked that, keeping Jh finite be- 
tween the e g electrons, but decoupling them from the ti g 
spins no longer gives rise to the FFF phase at low values 
of the JT splitting (cf. the phase diagrams in Figs. 6 and 
7). Instead at e ~ a FFA phase is found in agreement 
with the results of Ref. j|] for a model, which only consid- 
ered e g electrons and no JT splitting. This indicates that 
the determination of the stable phase (at least) at small 
values of the JT energy must take in due account the 
Hund coupling thereby including the t2 g levels. Secondly 
a quantitative determination of the stability region for the 
FFA phase and of the value of the magnetic couplings is 
subordinate to the consideration of the Jh term. 

Our work shares with Ref. (20j the generic result that 
JT distortions strongly affect the magnetic structure. Nev- 
ertheless it is worth pointing out some differences. In a 
certain respect our work is less ambitious in so far it does 
not attempt to determine the JT distortions, but it rather 



imposes them as external parameters of the calculation. 
Actually we do not believe that such deformations can be 
easily determined by microscopic models, which should 
incorporate complex effects such as long-range Coulomb 
interactions, cation and anion size and tilts of the MnO§ 
octaedra. On the other hand realistic deformations as ob- 
tained from experiments can easily be imposed and the 
consequent local electronic structure can be determined 
exactly: orbital ordering results essentially from coopera- 
tive Jahn- Teller deformations. 

Moreover, and quite importantly for a quantitative de- 
termination of the magnetic coupling and of the stability 
of the magnetic phases, we here also take into account the 
electronic Coulomb repulsion. This interaction is perforce 
larger than the JT interaction and contributes to its insu- 
lating behavior as well as to the numerical values of the 
exchange couplings. 

Finally we showed that using reasonable parameters 
the experimental values of the magnetic couplings can eas- 
ily be reproduced. Of course precise estimates depend on 
the knowledge of the various couplings entering the model, 
which are not always available neither from experiments 
nor from reliable first principle calculations. However cal- 
culating the magnetic couplings for various parameters 
and matching the numerical results with the experimen- 
tally obtained values provides useful connections between 
the involved parameters and set limits to the poorly known 
physical quantities. 
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